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AN ANALYSIS OF THE STABILITY OF AN AIRPLANE WITH FREE CONTROLS 
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SUMMARY 

An investigation is made oj the conditions essential to 
the stability of an airplane with free control surfaces. 
Calculations are based on typical airplane characteristics 
with certain factors varied to mver a range of eurrent 
designs. Stability charts are included to show the limiting 
values of the aerodynamic hinge moments and the toeight 
hinge moments of the control surfaces for various positions 
of the center of gravity of the airplane and for control 
systems with various moments of inertia. The ejects of 
reducing the chord and of diminatir^ the floating tendency 
of the surface, of changing the wing loading, and of 
decreasing the radius of gyration of the airplane are also 
indicated. An investigation has also been made of the 
nature of the motion of the airplane with controls free and 
of the modes of instabiliiy thai may occur. 

Stability with the controls free generally depends more 
critically on the design of the control system than on the 
stability characteristics of the airplane. In particular, 
too great a weight moment, combined with a high degree of 
aerodynamic balance, may cause undamped oscillations. 
Regardless of the weight moment, it appears difficfidt to 
secure stahUity when the aerodynamic balance exceeds 75 
percent of the hinge moment. 

INTRODUCTION 

During recent investigations by thB NACA of 
the flying qualities of several airplanes of different 
types, a tendency toward longitudinal instability was 
noted that involved pitching of the airplane reinforced 
by movements of the elevator. In other flight tests, 
lateral instability accompanied by oscillations of the 
rudder or of tlie ailerons has been noticed. The 
oscillations observed were rapid enough to be influenced 
by the inertia of the control surfaces but were not 
believed to be sufficiently rapid to involve the elasticity 
of the structure. The problem is thus concerned with 
motions intermediate between flutter and movements 
of the airplane as a rigid body. 

It was thought that a theoretical analysis of the 
stability of an airplane with the controls free might 
shed some light on the cause of these undesirable 
motions and might indicate how they could be avoided 
in design. Of the previous publications on the subject, 



the most detailed is that of E. Bartsck on lateral 
motions of an airplane with free rudder and ailerons 
(reference 1). In order to make specific reoommenda-' 
tions applicable to modem design, a study of stability 
more complete and detailed than any available was 
undertaken. Calculations were made covering both 
longitudinal and lateral motions and the elevator-frec, 
the rudder-free, and the aileron-free conditions. The 
computations were based on a set of typical airplane 
characteristics, except for parameters introduced to 
cover such variations in control-surface design as seem 
most likely to afiect stability. The results that might 
be expected under corresponding conditions in airplanes 
with different over-aJl mass characteristics have also 
been indicated. 

SYMBOLS AND COEFFICIENTS 

The following symbols are used in addition to those 
defined in the report covers. (See figs. 1 and 2.) 

The subscript c refers to a control-surface character- 
istic and is replaced in the various sections of the 
report by e for elevator, r for rudder, and a for ailerons; 
the subscript s refers to the control stick or wheel 
mechanism. 

A length equal to one-half the mean wing chord is 
used as the fundamental unit of length in order to obtain 
the results in a form applicable to geometrically similar 
airplanes of any size or loading. Conversion to this 
system is made by dividing all lengths measured in 
ordinary units by the length of the half-wing chord. 
Quantities entering into nondrmensional expressions do 
not, of coiuse, require such conversion. 

Z7o steady-flight speed 
j3 angle of sideslip 

ke radius of gyration of control mechanism about 

control-surface hinge axis 
Xc moment arm of center of gravity of control 

system about hinge axis, positive when center 

of gravity is behind hinge 
y distance from center of gravity of aileron to plane 

of symmetry 
A aspect ratio 
I tail length of airplane 
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projection on X axis of distance between 
center of gravity of airplane and its aero- 
dynamic center (with controls fixed) 



n =— ^airplane density ratio 



He= — ^control-surface density ratio 

{ aileron weight-moment parameter. (See 

equation (14).) 
H conti'ol-surface hinge moment . . 

1) control gearing ratio 

6=U4 distance along flight path 

. • . dw d5 . 
w, S, etc.=^.^) j^j etc. 

fl, S, etc.= ^> etc. 

D =^ differential operator 

C C C etc 

STABILITY WITH ELEVATOR FREE 

Pitching motions sufficiently rapid to be affected 
by the inertia of the elevator control probably will not 
involve sensible changes in the forward speed of the 
airplane. Accelerations of the airplane along the flight 
path will therefore be neglected. The rapidity of the 
oscillations makes it advisable, on the other hand, to 
include certain aerodynamic effects not retained, in. the 
equations of motion in their usual form. In addition 
to the moments developed in response to the displace- 




10, Z 

FiGCBE 1.— Notation for longitadinal mottoiis. 

ments of the tail surfaces, moments due to angular 
velocities of these surfaces are also considered. Thus, 
the pitching moment due to angular velocity of the 
elevator about its hinge bM/dh, the pitching moment 



due to the aerodynamic inertia of the surfaces bAI/bw, 
and the aerodynamic damping of the elevator bH/db will 
be included in the present analysis. Secondary factors 
entering into the equations, such as the vortical acceler- 




FiQUBE 2.— Notation for Istersl motions. 



ation at the center of gravity due to the lift of the 
horizontal tail, are ni^lected. The equations of motion 
take the following form: 



miw-Ui9)-w^ 

■ ,•• -bM bM .bM 

mkTd—6 — -—w^ w-x-r- 

i) 0 bw bv) 

bi dS 



,bM i^M 



— 0 

=0 

=0 



(1) 



If the following substitutions are made 



m 
2 



„ d l_d 

w=aUo 

1=' 



and the stability derivatives are replaced by the 
equivalent coefficients, equations (1) are reduced to the 
following nondimensional form : 
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=0 



=0 



=0 



(2) 



It is to be noted that, with the exception of the inde- 
pendent variable s, equations (2) involve no quantities 
dependent on the steady-flight speed. Motions plotted 
against tie distance s are therefore applicable to any 
initial flight condition within the unatalled range. 

The equations of motion are based on the assumption 
of a constant forward velocity and therefore do not 
show the possibility of a phugoid oscillation of in- 
creasing amplitude or tiie possibility of a ceirtain type 
6i slow divergence from the steady-flight attitude. 
Experience has shown that the unstable pKugoid motion 
is not likely to cause trouble under ordinary operafiDg 
conditions because its period is of the order of 2500 
chord lengths; the oscillations of interest in control-free 
stability have periods of the order of 50 chord lengths. 
The slow divergence corresponds simply to a loss of 
static stability when the control is free. Static stability 
with the elevator free is assured if the following con- 
dition is fulfilled: 

CkfCmJ^CmtCiig ■ (3) 

In this paper the divergence treated is of a more rapid 
type. 

Because the pitching that enters into the equations of 
motion was expected to be quite rapid, it was thought 
that the lag in the effect of the wing wake at the tail 
would be an important factor. Under steady condi- 
tions, the wing wake diminishes the relative angle of 
attack at the taU to about one-half. After a sudden 
change of angle, however, the tail will at first receive a 
strong upwash due to vortices shed by the wing in 
consequence of its additional circulation. The result 
is a rather complex transient variation of the vertical 
velocity. This variation affects both the lift of the 
horizontal tail and the floating moment of the elevator. 

The possible effect of the transient-flow phenomenon 
at the tail was estimated by making several calculations 
in which a simple fixed lag in the action of the downwash 
was assumed, expressed by setting ^ 

<^".-= 0^^+ O^^S^ -a^-) (4) 

0>,a=0.,„u^l-e^-"^ (5) 

A compai-ison of the resulting motions witii corre- 
spondiog results obtaiued when tiie lag function was 

1 For the use of the operator e-"> to show the effect of lafr, see leference 2, page !6. 
Subsequeat Investigation (refeience 3) has shonn the ramplez transient effect to Iw 
more nearly approximated hj the operator aae-C'^ 



entirely omitted showed that the lag, although having 
a noticeable effect on certain stable modes of oscillation, 
caused only a small change in the slower type of oscilla- 
tion jn which instability occurs first. Revision of the 
computations to include a more acciu^ate representation 
of the lag was therefore considered not worth while and 
all calculations were allowed to stand with egfT"* as 
the kg operator. In order to combine this operator 
with other tenns of the equations, the expression was 
expanded loto a power series in D, 

The stability of the motions is indicated by the nature 
of the roots of the characteristic equation, which is 
obtained from equations (2) by settiag the determinant 
formed from tlie coefficients equal to zero. If D6, 
rather than $, is considered one of the variables, this 
equation is 



—It 



—flfXe—Oitjj, 



0 



=0 



The equation is thus a quartic, and tenns introduced 
by the expansion of Gn„iP) and C^^{D) that would 
increase its degree were discarded because the roots are 
always small and higher powers are negligible in value. 

The roots of the stability equation were found for 
several typical cases. Apparently, in the usual case the 
motion is oscillatory and of two fairly distinct modes. 
One of the modes of oscillation, although more rapid 
than the modes encountered with the controls fixed, is 
nevertheless slow enough (with a period of the order of 
60 chord lengths) to involve coupling and reinforcing 
movements of the airplane. The damping is conse- 
quently light, and instability will occur first in this slower 
mode. It is undoubtedly this mode that has been ob- 
served in flight in the cases mentioned in the introduc- 
tion. The second mode is much more rapid but heavily 
damped. The short period (about 15 chord lengths) 
suggests that the motion is essentially limited to a 
flapping of the elevator and may become unstable only 
as flutter involving elastic deformations of the structure. 

It was expected that variations in the aerodynamic 
hinge-moment slope Cjj, the mass-moment coefficient 
tlie moment-of-iaertia coefficient iijc^, and the 
static stability coefficient Cn^ would be most important 
from the designer's point of view. These quantities 
were therefore retained in the equations as parameters, 
and numerical values were substituted for the remaining 
quantities. Limiting conditions for stability are tiien 
in the form of relations connecting the four variables. 

Of the conditions for stability, only two were found 
to be effective within the practical range of the param- 
eters. "A boundary beyond which straight divergence 
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occurs is obtained by setting the constant term of the 
stability equation 

equal to zero. This expression is independent of the 
elevator moment of inertia. The second boundary is 
the limit for oscillatory stability and is obtained by 
applying Routh's discriminant to the stability quartic. 
This boundary was found to shift only a negligible 
amount with a large change in the static-stability co- 
eflScient Cm„ and was therefore considered to be inde- 
pendent of Om^. Partial elimination of the parameters 
in tlus way made possible the presentation of the results 
in a simplified form. 

The computations of figure .3 were based on the char- 
acteristics of generally used types of balance and on a set 
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Mass -moment parameter. Mm ^ 

FiGVBB 3.— Elevator-free stability legloiis. (iO-petcent-<!hord elevator; ii, 4S; kr, I.TB. 

of typical airplane characteristics. (See table I of the 
appendix.) The effect of variation of the moment of 
inertia was subsequentiy investigated. Figure 4 covers 
the case of an airplane with the radius of gyration re- 
duced to make the moment of inertia half the average 
assumed for figiire 3. An investigation was also made of 
the stability of a more heavily loaded airplane by dou- 
bling the density factor n and comparing (fig. 5) a repre- 
sentative stability boundary (in terms of nJc^ and n^,) 
with the corresponding curve for the conditions of 
figure 3. The particular variations chosen were con- 
sidered representative of the trends in modem airplane 
design. 

Of the over-all characteristics of an airplane, the 
radius of gyration seems most likely to affect its 
stability. The results show that an airplane with a 
small radius of gyration will not permit so wide a range 
of the elevator design parameters as will the assumed 
average airplane. Its greater responsiveness to elevator 
deflection will c-ause it to reinforce more readily the 



movements of the elevator leading to oscillatory in- 
stability. The boundary for divergence (equation (7)) 
is independent of variations in the radius of gyration. 
As shown in figure 5, the relative density or loading 

r 
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Mass-momenf parameter, 

FiQUKE 4.— JSlevator-frce stabfltty. JtegloDs for reduced airplane moment of Inertia . 
fO-percect-chord elevator;jt, 46; kr, US. 

of the airplane is not a critical factor, which may be 
attributed to the fact that the normal relation between 
the lift-curve slope and the loading is such that the 
airplane is effectively constrained against relative 
motions normal to the wing surface. Differences in the 
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FiGUXS 6.— Elevator-f^ stability. ESect of increasing density factor. Ct,, ; 
I...., 0.12; kr, 1.79. 

degree of this constraint, as caused by ordinary varia- 
tions ia either or /t, are unimportant. 

In general, itTnay be concluded that the design of 
the elevator itself is of critical importance in obtaining 
control-free stability. A large mass moment or moment 
of inertia of the control surface is seen to be unfavorable 



AN ANALYSIS OF THE STABIUTY OP 



AN AIRPLANE "WITH FHEE CONTROLS 



113 



to stability. Of primary concern, however, is the ad- 
verse effect of aerodynamic balance, especially because' 
it is found necessary to resort to a high d^ee of balance 
with many modem airplanes. 
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Aerodynamic balance area, percent 

FiouBE S.— Typical TSriatlon of Ct, and with aetodynamlc balance area fiii small 
deSecttons (trom reference 4). 

Figure 6 (taken from data of reference 4) shows typical 
hinge-moment-coefBcient curves for a control flap having 
the. inset-hinge type of balance used in most modem 
control systems. In these experiments, tJie hinge mo- 
ment due to a unit change in the angle of attack 
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Mass- moment parameter, ti^Sg 

FiacRB 7.— Elevator-free stabUIty. Begions for nonfloatlDE elevator. SO-peicent- 
chord elevator; a 41; kr, l.TB. 

remained practically constant as the balance area was 
increased^ This form of balance thus would not provide 
compensation for the floatiag moment G^^ in the same 
proportion as for the restoring moment C*, and, as the 
degree of balance was increased, the equilibrium floating 



angles would become increasingly large, so that there 
would be greater danger of static instability with controls 
free, as shown in equation (3). The same considerations 
apply to the bfllancing tab. 

On the other hand, it should be possible to com- 
pensate for the floating pressure in the same or, perhaps, 
in a greater proportion than the proportion of reduc- 
tion of the restoring moment. Thus, with a horn type 
of balance, for example, the equilibrium floating angles 
may be held constant or may even be reduced, which 
results in greater static stability.* A comparison of 
figure 3 (Ca = — 0.24) with figure 7 (C* =0) shows 

"tetl ^ ^ "/oil 

that decreasing the floating moment also decreases the 
likelihood of rapid divergence. The boundary for 
oscillatory stability is hardly influenced by this factor. 

The computations for either type of balance apply 
to an elevator operated by a servo tab, provided that 




-4 0 ' 8 la 16 
Mass -moment parameter, ft,\ 

FionaE 8.— ElevatOF-free stability. Regions fox redaced elevator cbord. i!S-percent- 
cbord elevator; n, 45; ir> I.7B. 

the tab remains fixed relative to the elevator during the 
oscillations. Thus, as far as stability is concerned, 
servo operation with controls fixed corresponds to the 
ordinary control-free condition. The stability with 
both elevator and servo tab free is not covered in the 
present study. 

Future designs wiU probably show a trend toward 
narrower control surfaces, whether balanced or not, 
because the basic hinge moments can be markedly 
reduced mth a small loss of effectiveness. If the chord 
of the elevator is reduced from 50 percent to 25 percent 
of that of the horizontal tail surface, its effectiveness is 
reduced by only 30 percent; whereas, tlie basic hinge 
moment is divided by 4- 

Figure 8 shows the regions of stabiUty mth a reduced 



> Anotber advantage oC the horn type at balanca Is Ibat the hinge gap may be 
sealed. The anbject of horn halancea ts discussed (Urtber by HempbiU in refarence 5. 



114- 



BBPORT NO. 709 — NATIONAIj ADVISORY COMWITTBE FOB AEEONAtJTlCS 



elevator chord. The differences between, these r^ons 
and those of figure 3 are principally due to changes in 
the coefficients bCJba, dC^/bDS, and (iOJbDS. As 
previously noted, the control moment (proportional to 
5(7m/£>5) is reduced hy only 30 percent. In the inter- 
pretation of this figure, it should be borne in mind that 
the ratios x, and k, would naturally be smaller for the 
narrower elevator. If account were taken of this scale 
factor, the region of stability would appear much wider 
than the region for the SO-percent elevator. 

An effective method of obtaining greater stabilitj m 
the control-surface motions is tiie introduction of 
additional damping into the syston. If the responsive- 
ness of the control surface is reduced, a considerably 
larger degree of aerodynamic balance may be used 
(fig. 9). The permissible mass unbalance is also 
increased, although to a lesser extent. The results 
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9.— Blevator4ree stabfUt; regions. Effect of additional damping In the 
csontrol system. EO-percentHdioid elevstor; ^ i&\ kr, 1.79, 



shown are for a comparatively smaU amount of damping 
(AC»^,= — 1.0), which corresponds, for a rate of 

deflection of 20° per second, to the force required of 
the pilot to maintain 1° of elevator deflection, 

STABILITY WITH RUDDER FREE 

Because the lateral motions involve two controls and 
five degrees of freedom, the analysis is more complex 
than for the longitudinal motion, which has one control 
and three degrees of freedom. Fortunately, the 
rudder and the ailerons exert, their principal influences 
on different modes and only a slight loss in accuracy is 
incurred if each mode is treated separately. 

Oscillation of the rudder control wiU. be primarily 
influenced by coupling with the yawing oscillations of 
the airplane. The small rolling oscillations simul- 
taneously induced will generate neither very strong 
yawing moments nor very strong rudder hinge moments ; 



hence, the rolling degree of freedom will be neglected 
in the.._examination of rudder-free stability. This 
assumption, which has been checked quantitatively by 
Bartsch (reference 1), reduces the simultaneous equa- 
tions of motion to the following form: 



mUo(r+P)-^ 



=0 



m „m ^bN ;dJV 



S- 



bp bS 



-S- 



bS 



bH bS 



_lbH_:bH 
b« bf 



=0 



(8) 



If substitutions corresponding to those introduced in 
the elevator calculations are made, the equations are 
reduced to the following nondimensional form: 



(9) 



The stability equation is 



ft 



'»/» 



nXD+0,p iirik*+x^D 

— (C»j,#+A*,Sr) 



-(C.wO+C.,) 
ihkr'D" " 



=0 



Thia^quation is closely analogous t0-eqtration (6) for 

longitudinal motion. The corresponding coefficients 
have similar values with the exception of (7y^, which is 
much smaller than the corresponding term Cl^ because 
the normal force that is developed by the wing in 
pitching is absent in the lateral motions. 

Thexoots of the stability equation again indicate two 
modes of motion. Thus, in a typical case, the roots 
are — 0.008±0,D35i and — 0.25±0.28i.. In this in- 
stance, the modes are both oscillatory. The first 
pair of roots indicates a lightly damped oscillation oL 
such a frequency (period— 90 chord lengths) as to involve 
sensible coupling between the yawing of the aii-plane 
and the swinging of the rudder. The second mode is of 
much higher frequency and undoubtedly represents the 
natural oscillation of iJie rudder with the airplane acting 
as a practically rigid support. When the restoring 
moment C*, is reduced, the second mode becomes 
aperiodic and eventually divergent as the motion bc- 
comes'less rapid. Oscillatory instability appears first 
in the slower mode, as in the case of the elevator. 



Page 115, figures 10 to I3: All values of O^g given in figures 10 
through. 13 should he reduced hy dividing by 6; that" is, 
the values on the curves shoijld read Cng = 102, ,OSh, 
and ,026 instead of .6l2, and .155. 

Page 116: In the second equation of the two "bracketed as (l2) , the sign 

of the third term should he changed from minus to plus; thus 
the equation should readl 



, 3 d^6 as , dT3/ — , 



Pages ll6- and lljl Figure lU, containing the aileron-free stability 

"boundaries, has heen found to he incorrect and 
shotild not "be used. The paragraph referring to 
this figure, "beginning at the "bottom of page II6 
and continuing through equation (15) on page 117» 
is therefore also ia error. The correct "boundaries 
for oscillatory stability with ailerons free are 
given in a new figure 1^ which has been inserted 
in this copy of the Annual Report. Ihe boundary 
for divergence is that given by equation (16) . 
Statement 3 under "Concluding Eemarks," page 117» 
no longer applies to ailerons. 




Aerodynamic hings-mornent parameter, Cj,. 
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Cnp. The greater the value of 0„p, the less is the 
allowable mass moment with any given amount of 
aerodynamic balance. On the other hand, if the 
mass moment is small enoiigh to insure damping of 
the oscillations, a larger value of C7„^ will increase the 
aerodjmamic balance that may be introduced without 
causing divergence. The radius of gyration of the 
airplane is of considerable importance, shifting the 
boundary for oscillatory stability so as very nearly to 
double the stable region when the moment of inertia is 
doubled. (See%. 10.). 

Instability with the rudder free is likely to occur in 
the form of a divergence. The criterion is practically 
the same as the condition for weathercock stability 
with the control free, which is . 

c*A^o (11) 

This critei'ion is independent of the moments of inertia 
of the airplane and of the control surface. The greatest 
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gain in this margin of stability is obtained by increasiog 
and reducing the floating tendency of the rudder. 
Eeducing the chord of the rudder decreases both Oh^ 
and Cnt and hence considerably widens the margin of 
stability. (See fig. 12.) Complete elimination of C*^ 
by the use of a nonfloating type of balance eliminates the 
likelihood of divergence within the normal range of 
weight distribution (fig. 13). Such a gain, however, 
would be achieved only by sacrificing some margin of 
oscillatory stability. 

The lateral oscillations of an airplane witli controls 
fixed are known to be influenced by coupling between 
tlie rolling and tlie yawing motions. These oscillations 
tend to become undamped when the weathercock 
stability approaches zero. Rrewng the rudder 
control diminishes C'ajj and may thus lead to this type 
of oscillatory instability. The condition for zero 
weathercock stability is approximately that for straight 



divergence in the control-free condition (equation (11); 
and the boundaries for divergence (figs. 10 to 13) can 
therefore be interpreted also as boundaries for stability 
of the slow lateral oscillation, when the airplane is free 
to roU, The absence in the criterion of terms involving 
the mass of the rudder may be explained by the fact 
that, as the limit of stability is approached, Uie osciQa- 
tion becomes very slow and tlie ya\¥ing component lends 
to disappear. 

STABILITY WITH AILERONS FREE 

The stability of an airplane with the ailerons free is 
examined by including in the equations the interaction 
between, rolling motions of the airplane and movements 
of the ailerons. Small simultaneous yawing and side- 
slipping motions will also occur but, because their 
reactions on the rolling and the hinge moments are 
small, they may be neglected. The resulting equations 
are: 



'""dt-Pb^-^bs-^-^. 



=0 



d6 



or, in nondimensional form. 



■bH 



(12) 



=0 



(13) 



In the aileron control, a part of the mechanism 
normally rotates about an axis at right angles to the 
aileron hinge axis. As a result, a part of the aileron 
torque produced by angular acceleration of the air- 
plane is proportional to the product of inertia of tlie 
aileron itself and another part is proportional to the 
moment of inertia of the control stick or wheoh . Both 
quantities are included in the parameter ^ Thus: • 



S^ea 



(14) 



where the subscript s refei-s to the control stick and ij 
is a constant inserted to take account of any difference 
of gearing between the control and the aileron. 

The calculations, based on values of the derivatives 
given in the appendix, cover two aUerojn widths: 15 
percent and 30 percent of the wing chord. Jn both 
cases, the ailerons were assumed to cover 50 percent 
of the wing semispan. 

Variations in the floating tendency of the aileron 
were also considered but were foimd to have little effect 
on the stability. The results given may therefore be 
applied to any of the existing types of balance that 
give smooth hinge-moment curves." _^ 
^he boundaries for stability in the two possible 
modes are presented in figure 14^ Instability appears 

* Certain ailerons oC the Frise type tbat show a reversal o( the hinge-moment slope 
may develop unoontrollable oscillations of fixed amplitude. This condition Is dls- 
cossed tn reference 6. 
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in tke form, of increasing oscillations and is most likely 
to occur when the ailerons have a high degree of aero- 
dynamic balance and too great a mass moment or 
moment of inertia. As the moment of inertia decreases, 
! the boundary for undamped oscillations approaches 
' that for straight divergence (which is itself independent 
of the moment of inertia); and it becomes possible for 
both types of instability to eixist simultaneously. The 
condition to prevent straight divergence, 

C„^-C,^^C;^,-G,,^O,,,-:^C,,^^'>0 (15) 

i 

,. is, however, always satisfied when there is damping of 

V. the- oscillations. 

Straight divergence, as encountered in the elevator- 
free and the rudder-free conditions, is indicated when 
the constant term of the stability equation is n^ative, 
that is, when 

Ct„^C,,^Ct^Ci, (16) 

This condition is not likely to occur unless the aero- 
dynamic balance is nearly complete. 

CONCLUDING REMARKS 

Experience has shown that, before the actual limit 
of stability is reached, tlie airplane undeirgoes oscilla- 
tions which, although damped, are still persistent 
enough to be undesirable. The boundaries given in 
the stability charts are therefore of value chiefly as 
indications of the effect of certain design factors; they 
are useful quantitatively only as outside limits, not to 
be approached too dosely. Further experiments will 
be necessary to determine the mai^^ of stability re- 
quired for smooth operation in gusty air. 

On the other hand, the charts are to a certain extent 
conservative because they do not take into account the 
possibility of friction of the control system, a factor 
that would widen the margin of stability. 

The indications of the present study may be sum- 
marized as follows: 

1. There is a limit to the effectiveness of the aero- 
dynamic balance that may be safely employed with 
any conventional control system. In most cases, it 
appears difficult to secure stability with the hinge 
moment reduced to less than 25 percent of its value 
for the imbalanced surface. 



2. Reduction of the floating moment Ciuc, if it can 
be brought about independently of a reduction of the 
aerodynamic balance, causes a shift of the boundary 
for divergence. (Cf. figs. 3 and 7 and figs. 12 and 13.) 
The effect is particularly noticeable in the case of the 
rudder, where the likelihood of this form of instability 
is materially decreased. 

3. Within the usual range of .characteristics, the 
elevator a nd -Uia JaflBrDg-contrQjg^^ more susceptible 
to oscillatory instability than to the rapid form of 
divergent instability. The stability with -eit-her-of V^<'j 
-these- control! free may be improved by (a) using a less 
effective aerodynamic balance, (b) decreasing the mass 
moment and the moment of inertia of the control, or 

(e) using a control surface of narrow chord. 

4. Divergence is a more likely form of instability for 
the rudder control (figs. 10 and 11) and may be avoided 
by reducing the effectiveness of the aerodynamic bal- 
ance or, as has been suggested, by using a balance that 
reduces the floating tendency of the rudder, although a 
highly effective balance of a type that reduces the 
floating tendency may result in oscillatory instability. 

6. The oscillatory stability of the elevator-free 
system is but little affected by the restoring moment of 
the airplane in pitch {G„J. In the case of the yawing 
motions, however, the existence of a strong restoring 
moment (Ca^) increases the likelihood of oscillatory 
instability. (See, for example, fig. 10.) 

6. In all cases, an increase in the relative radius of 
gyration of the airplane results in an increased range of 
stability (cf. fip. 3 and 4; see also fig. 10), but changes 
in weight without corresponding changes in the rotary 
inertia have little effect. (See fig. 5.) 

7. The use of a narrow control surface is recom- 
mended as a means of increasing the control-free 
stability as well as from other considerations. The 
marked effect of reducing the chord is sho^vn by a 
comparison of figure 8 with figure 3 and of figure 12 
with figure 10. 



Langlet Meaiosial Aeronautical Labobatoby, 
National Advisobt Committee fob Aeeonautics, 
Lantglet Field, Va., August 16, 1940. 



APPENDIX 



STABIUTV DEKIVATIVES 

The geometric and the aerodynamic characteristics 
used in the stability calculations are given in tables I 
toIV. ' . . - 

Tablh L— general airplane CHARACTERISTICS 



Tall length j^yj 

Wing chord " " ' " - - r ..- . 

Horizontal-tan area iks 



Wing area 
Vertlcal-tall area 




06 


Wing area 
Horliontal-tall chord 




,..„ .i!0 


Wing chord 






Vertlcal-tall chord 




3S3 


Wing chord 
Aileron chord 




IS 


Wing chord 
Aileron span 




60 


Wlog semispan 










. 1L0(1 



Table II.— ELEVATOR-TREE STABILITY • 
COEFFICIENTS . 





SO-percent-ohotd 
etevBtor 


26-percent-chotd 
elevator 


^ - 


-is 

-9.26 
4.8(ai8«-j;.,.,) 
-1,121 
-t4l» 
-.060 
-.57 
• -1.00 
-I.8B 
-a24aDdO 


4.8 
-Has 

4.8(0.185-*....) 

-1.121 
-I.45D 
-.672 
-.28 
-.50 
-.408 
-.075 


ciy 

*«*»« 




"'Dn 






X"" 




OlD, 


"un 





Table III.— RUDDER-FREE STABILITY COEFFICIENTS 





SO-perCent- 

chord rudder 


25-perceiit- 

chord rudder 


p -'- 


-a 274 


-0.274 




-^682 


-.682 


^- 


-.0567 


-.0897 


X"ci 


—.0460 


-.0176 







—.276 





.123 


.075 




-1.0O 


-.50 



Tablh IV.— AILERON-FREE STABILITY COEFFICIENTS 





15-pereent- 
chord afleton 


30-peccent- 
chord alleion 




C"! 




— 

-n.lW 

-.184 
—.680 


-2.94 
—.352 
-.296 
-.322 

-1.S87 



The aerodynamic coeflBcients arc, in most cases, 
based on experimental results; theoretical values are 
used only where such results were not established. 
Discussions of the more commonly used derivatives will 
be iound in references 7, 8, and 9. Several of the 
unfamiliar coefficients are developed in the following 
paragraphs. 

Damping in pitching CW^^. — The piincipal com- 
ponent of damping in pitching, furnished by the horizon- 
tal t^, is 



In addition, the pitchmg motion introduces a relative 
camber of the. wing section, giving rise to a moment . 
coefficient 

4 2C7o 

Thus 

Then, since 

Uo e 

- — - e„„,— [2|(L)K,+i] 

Htching-moment slope Cm^. — ^The pitching-moment. 
slope is given by 

where Xa.c. and Za.c. are the distances of the aero- 
dynamic center of the complete airplane beliind and 
below the center of gravity. The location of the aero- 
dynajnic center for the airplane as a whole is estimated 
by taHng the centroid of the aerodynamic centers of the 
various components. Thus, if terms in Za.e. (which is 
usually small) are neglected. 
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Fitchmg-moment coefficient due to vertical accelera- 
tion Cmj„. — The pitching-moment coefficient due to 
vertical acceleratioa arises from the aerodynamic 
inertia of the wing and the horizontal tail surfaces in 
motion normal to their chords. The force on each 
surface is equal to the reaction of a body of air described 
by rotating the surface about its midchord line. Thus, 
the pitching moment 



M 



63:50+ TP 



and 




where is the distance from the 50-percent-chord 
point of the wing to the center of gravity of the airplane 
and siiOf is the same distance measured from the cor- 
responding point of the horizontal tail. 
Pitching moment due to elevator deflection Cg,,. — 

The pitching moment due to elevator deflection is 
given by the formula 

St I 



where a, is the angle of zero lift of tJie devator. 

Theoretical and experimental values for ikejbd, for 
flaps witii sealed hinges, are given in figure 15. The co- 
efficients used in these calculations ((7,, aad as veil 
&s Om^) were, however, based on the experimentally de- 
termined changes of lift produced by a flap with open 
gaps at the hinges. The effect of a small gap is to reduce 
the effectiveness of the flap by about 30 percent. At a 
large deflection, the flap with inset-hinge balance shows 
a stUl greater loss because of the protruding balance 
portion. 

Fitching-moment ooeffldent due to angular velocity 

of the elevator C'„^j. — The pitching-moment coeffi- 
cient due to angular velocity of the elevator is 



where 



dat 



The parameters ^Cr,J(>D5 and ^a,/bDd may be found as 
functions of the chordwise position of the hinge from 
figure 15. The figure is based on the theoretical treat- 
ment of Theodorsen (reference 7), with the assumption 
of long oscillations ^eater than 20 cbord lengths). It 
must be remembered that DS involves the distance 
traveled by the airplane measured in terms of its half- 
wing chord, and the quantities given must be multiplied 
by the ratio Ctfca to convert them to half-wing chord 
lengths. 

Damping moment of elevator Cj^j. — The hinge mo- 
ment due to angular velocity is treated theoretdcaUy 
by Theodorsen in reference 7. Figure 16, derived from 



the theory, gives the component parameters of the 
damping moment as functions of the chordwise position 
of the hinge. The same considerations are effective 
here as in the apphcation of figure 15. 

Hinge moment due to pitching C7j^. — Positive pitch- 




Fiohbb is.— Fanmetets toe detennlnlng the effects ot aognlar Telodtj and deflection 

of flaps on thB lift. 

ing motion causes an increase in the angle of attack of 
the tail surface equal to Is/Uq. The resulting hinge- 
moment coefficient is 

IB, 



and, since 



TJm c 
c "urn 

Inasmuch as rotation of the airplane about its center of 
gravity does not appreciably change the lift of the 
wing, the downwash correction may be neglected. 
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Hinge moment due to change in angle of attack 
C^^. — In accordance with the indications of testa by 
Goefct and Reeder (reference 4), the aerodynamic float- 
ing moment C^^^ was assumed independent of the 

degree of aei'odynamic balance of the control surface 




FiGUBE 16.— Parameten for detennlnlng the dampliig moments of lisps. 



(fig. 6). . The assumption is valid for the inset-hinge 
type of balanced flap shown in this figure. The floating 
moment will vary with the type of balance, however, 
as discussed in the text, and additional computations 
were tlierefore made in which (7> ' was assumed equal 

to zero. 



The lateral-stability derivatives Cr^, 
50. ^ _ i>0. 



a 



110, 



^"Di- 



re 



pe 



are discussed in refer- 



2Uo 2Uo 

ences 8 and 9. The other coefficients for the lateral 
motions are derived in a manner closely analogous to 
the derivation of the corresponding longitudinal 
coefficients. 

The values of the mass moment and the moment-of- 
inertia coefficients of several representative elevator- 
control systems were determined experimentally in 
order .to find the magnitude and the range to bo expected 
in practice. The experiments were made by attaching 
a spring of known stiffness to the control column, oscil- 
lating tlie system, and recording the frequency and the 
damping. The mass moment was measured directly 
with a spring scale. The interpreted results are given 
b the following table: 



•Airplane 








2.70 


1.3JS 


Coitlss P-8« 


0 


.SO 




0.6 


4.5 


Paircbild 22 


1.0 


LI 
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